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Abstract 

We present analytic solutions for steady flow of the Johnson-Segalman (JS) model 
with a diffusion term in various geometries and under controlled strain rate condi- 
tions, using matched asymptotic expansions. The diffusion term represents a singu- 
lar perturbation that lifts the continuous degeneracy of stable, banded, steady states 
present in the absence of diffusion. We show that the stable steady flow solutions in 
Poiseuille and cylindrical Couette geometries always have two bands. For Couette 
flow and small curvature, two different banded solutions are possible, differing by 
the spatial sequence of the two bands. 



1 Introduction 



Many experimental results confirm the possibility of shear banding, i.e. the 
separation of bands of different shear rates and apparent viscosities in the flow 
of various systems of wormlike micelles (aqueous solutions of surfactants and 
salt [1-5] or organic solvent solutions of metallic complexes [6]), or lyotropic 
liquid crystals [7]. This type of behaviour can be explained as being the result 
of a constitutive instability, suitably described by non-monotonic flow curves 
such as those arising from the Johnson-Segalman (JS) [8] or Doi-Edwards 
models [9]. It is possible that similar constitutive instabilities are responsible 
for spurt and extrudate distortions of polymer melts [10]. 
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The steady banded flow solutions of the JS model were previously studied in 
planar [11,12], Poiseuille [13-17], and cylindrical Couette geometries [18]. In all 
these geometries banded flow solutions have continuous degeneracy [16,17,19], 
resulting from the indeterminacy of the positions of the interfaces and of the 
number of bands. The stress value can be kinetically selected [20,12,18,16,17], 
and it may differ from the top jump prediction of Refs. [13,14], although the 
selected value depends on the flow history and on the imposed shear rate 
[19]. This is in contrast to experiments on wormlike micelles that show a well 
defined stress plateau in the flow curves: the total stress is history independent 
and do not change with the shear rate. Furthermore, the evolution of the stress 
during transient flow may be rather long, suggesting the slow migration of the 
interface between bands to an equilibrium position [2]. 

Recently [19,21,22] we showed that, by supplementing the JS model with a 
stress diffusion term, one can lift the degeneracy of the steady flow and also 
account for the slow migration of the interface. The origin of this diffusion 
term can be justified in terms of the Brownian movement of polymer chains 
in an inhomogeneous stress field and can be deduced from the Fokker-Planck 
equation for a system of dumbbells [23]. There is no calculation, at present, 
for such a term in the wormlike micelle system. Similar results have been 
obtained by Refs. [24,25] by introducing phenomenological diffusion terms in 
toy constitutive models. Although non-local terms are not important in ho- 
mogeneous fiow, in banded fiow a realistic description of interfaces as strongly 
inhomogeneous regions requires such terms. This necessity has been noted by 
many authors [23,26-28]. 

Recent experiments [29] on wormlike micelles in pipe (Poiseuille) and cylindri- 
cal Couette flow geometries, or in cone and plate [30] using NMR microscopy 
found shear rate proflles with two or three bands. Birefringence measurements 
in Couette geometry [31,32] suggest the existence of only two bands, although 
the relation between birefringence and shear rate is hard to establish because it 
depends on microstructure. Most experiments on shear banding are performed 
in geometries imposing inhomogeneous total stress (Couette, pipe flow, cone 
and plate). Typically a stress quantity {e.g. torque, or pressure drop, or shear 
stress) is measured as a function of a rate of flow quantity {e.g. mean strain 
rate, flux through a pipe, wall speed), and an engineering "flow curve" is con- 
structed. Under inhomogeneous flow conditions the flow curves are not simply 
the local constitutive relation of the fluid. It is therefore important to assess 
theoretically the question of how many bands can occur and which are the 
flow curves for a given constitutive model and in various geometries. 

We present elsewhere numerical computations of flow curves in Couette flow 
[19] for the JS-d model. Here we show how matched asymptotics techniques 
can be used to obtain analytic results on the number of bands, interface width, 
and flow curves in various geometries. The application of these techniques to 
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the JS-d model represents only an example, and the same method can be 
easily adapted to other constitutive models. 

The structure of the paper is the following. In Section 2 we introduce the flow 
equations in various geometries. The matched asymptotic solutions of these 
equations are presented in Section 3. In Section 4 we discuss how the inner layer 
solution controls the history independence and stability of the flow. We apply 
our results in Section 5 to predict the shapes of flow curves and the interface 
width. In Section 6 we discuss the results and suggest possible extensions. 



2 Flow equations in veirious geometries 

2.1 General equations 

The dynamics of the diffusive Johnson- Segalman fluid is described by two 
equations. The first is the momentum balance: 

p(at + vV)v = V-T, (2.1) 

where p is the fiuid density, v is the velocity field, and T = — p I + 277D + S is 
the total stress tensor where the pressure p is determined by incompressibility, 
S is the viscoelastic stress carried by the polymer strands, 1] is the "solvent 
viscosity" and D is the symmetric part of the velocity gradient. The second 
equation is a constitutive relation for the polymer stress E, which we take to 
be the the JS constitutive model [8] with an added diffusion term [21,19]: 

S = W^S + 2-D--E, (2.2) 
T r 

where is the "polymer" viscosity, r is the relaxation time, and T) is the 
diffusion coefficient. 

The time evolution of S is governed by the Gordon-Schowalter (GS) time 
derivative [33], 

S = (^t + V V) S - (riS - S O) - a (DS + ED) (2.3) 

where a is the "slip parameter" and ft is the antisymmetric part of the velocity 
gradient. For a — 1 there is no slip and the GS derivative becomes the upper 

convective derivative. For a = (total slip) the fluid stress is not transmitted 
to the polymer strands and these can only be oriented by the flow: the GS 
derivative becomes the corotational derivative. In this work we consider \a\ < 
1. 
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We study the case of a fixed average shear rate, which represents a global 
constraint on the velocity gradient. We discuss several shear geometries: pla- 
nar shear, slit and pipe Poiseuille flow, and cylindrical Couette flow between 
concentric cylinders. We consider parallel stream lines in all geometries and 
leave the possibility of secondary flow for further investigation. 

2.2 Planar shear and Poiseuille flow 

Parallel stream line flow corresponds to v = v{ii)ex, for planar shear between 
plates at positions y = 0,L and for Poiseuille flow through a rectangular slit 
of infinite length and walls at positions y = ±L. 

We introduce dimensionless variables, V — —Vj^\/1 — a^W — {^^Hyy — 



^I^'/ = S^^^'^ = Ty^ly/l^.i = t/r.r = y/L. V is the velocity at 



7/ = (the maximum absolute value), having opposite sign to the shear rate 
and stress. With these definitions V has the same sign as the reduced shear 
rate and stress. In terms of these quantities the momentum balance reads: 



The dimensionless pressure gradient / vanishes in planar shear, e = ^ is 

the retardation parameter (viscosity ratio) and a = £^ is the ratio of the 
Reynolds and the Weissenberg numbers. In shear-banding experiments {e.g. 
for wormlike micelles) a ~ 10~^ — 10~^, therefore in this case one can neglect 
inertia and replace Eq. (2.4) by its creeping flow limit: 



adiV = dfS + edlV - f. 



(2.4) 



(2.5) 



The solutions of Eq. (2.5) are: 



a{r) — const 
a{r) = fr 



(planar shear) 
(slit Poiseuille). 



(2.6a) 
(2.6b) 



where a :— e'j + S is the total stress. 



The constitutive equations read 




(2.7a) 
(2.7b) 
(2.7c) 



dfW = DdlW -W + ^S 
d{Z = Dd^Z - Z, 
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and the shear rate condition is: 



V 



i: 



'y{f)dr. 



(2.8) 



2.3 Cylindrical Couette flow 



In the Couette geometry with concentric cyhnders at radii Ri < R2 and 
assuming circular stream Unes we have v = v{r)e.Q. The momentum balance 
in cylindrical coordinates reads adiv — ^df[f'^{S + £7)], where 7 = f-^ 

is the shear rate and a = ^ has the same significance as in planar shear. 
Creeping flow corresponds to: 



67 + 5 



(2.9) 



The constitutive equations read 



dtS 



DA.S - S 



AD 



dfW = DAfW -W + 



diZ = DA^.Z 



(1 — a^)r^ 
AD 

7 



(1- 



s 



[1 — a'^)f'^ 



+ 



ADa 



'1 — a'^)f'^ 



W 



where — df + \df is the Laplacian. The shear rate condition is: 



"i+p . dr 
r 



(2.10a) 
(2.10b) 
(2.10c) 



(2.11) 



where p — ^'^j^^^ ■ 

All the above variables are dimensionless and have the same sign, with the 

l+g 
2 



following definitions: V 

S 



I l+ay \ 



— a 



W 



^rr)~5 Z 



Tre^Vl — a^, and f = r/Ri. T is the torque per unit length applied at the 
inner cyhnder, and has the same sign as V, the velocity of the inner cyhnder, 
being opposite to the shear rate and stress. With these definitions T and V 
have the same sign as the shear rate and stress. 
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2-4 Pipe flow 



Like for Couette flow we use cylindrical coordinates, with the z axis parallel to 
the stream lines, v = v{r)ez- The pipe is an infinite length cylinder of radius 
R. The momentum balance reads adfV = ^df[f{S + 67)] — 2/, where 7 = |f 
is the shear rate. Creeping flow corresponds to a total stress that is linear in 
f, like for the slit geometry: 

+ S ^ a{f) ^ ff. (2.12) 

The shear rate condition is: 

t>= [^^(r)df. (2.13) 
^0 



The constitutive equations are: 



dfS = DAi^S - S 



1 + 



D 



+ 7(1- W^) 

ADX 



dfZ = DAfZ -Z + 



DX 



dfX = DAfX + X 



1 + 



4D 



We have used the following rescahngs: V — — o?, W 



(2.14a) 

(2.14b) 
(2.14c) 
(2.14d) 



l-a 



2 "^^2 



l+g V \T -y- 



/ 



1 dp rVT^ 

2 dz fj,T 



(E,,-E,,)^, Z = (i=«E,, + i±^E,,)^, 5 = E^.^v^r^, 
, D = P;^, and f = r/i?. 



3 Matched asymptotic solution for the steady flow 



The flow is described by the system of nonlinear, parabolic partial differen- 
tial equations of the reaction-diffusion type, Eqs. (2.7), (2.10), or (2.14). The 
nonlinear reaction terms are due to the polymer stress relaxation. The non- 
affine deformation (slip) is essential for the nonlinearity, and all rescalings are 
possible for \a\ < 1. 

We consider small D, for which the diffusion terms represent a small pertur- 
bation to the steady flow equations. Nevertheless, this perturbation is singular 
and the solution can not be represented as a uniformly convergent power series 
in D. The usual technique applying to this situation is asymptotic matching 
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[34]. The solution is divided into an inner layer solution around the interface 
where the diffusion terms are important and outer layer solutions a distance 
from the interface farther than its width (that scales like D^^'^), where diffusion 
terms are exponentially small and can be neglected. 



3.1 Outer layer solution 



Neglecting terms of order D and imposing stationarity in any of either Eqs. (2.6) 
and (2.7), or Eqs. (2.9) and (2.10), or Eqs. (2.12) and (2.14), depending on 
geometry, we obtain the same equations: 

t a{r) - Sout . 
lout = (3- la) 

Sout = ^^^(1 - Wout) (3.1b) 

Wont = ^o. /^"^ ~ '^""^ (3.1c) 
Zout = 0. (3.1d) 

For pipe flow we may also use Xout — Wont that together with Eq. (3. Id) 
implies {^e,e)out = 0- 

The algebraic system of Eqs. (3.1) is parametrised by a, the value of the total 
stress at the position f in the gap. The shear rate 7, as a function of the 
total stress a represents the local constitutive curve shown in Fig. la. Banded 
flow is possible between the minimum {a = o'^ottomi^)) and the maximum 
{a — crtop(e)) of the flow curve, when for one value of the stress there are three 
possible solutions of Eqs. (3.1). Only two of those (low {7 (a), S~{a), Z~{a)} 
and high shear rate {7 {&) , {&) , {&)}) correspond to stable bands, the 
intermediate shear rate being unstable [13,14]. The dependence of dbottom and 
atop on e can be obtained from Eqs. (3.1), and is represented in Fig. lb. 



3.2 Inner layer solution 



In any of the Eqs. (2.6) and (2.7), or Eqs. (2.9) and (2.10), or Eqs. (2.12) and 
(2.14), we impose stationarity and change variables to: 




where f* is the position of the interface. 
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Fig. 1. a) Local constitutive curves in reduced variables, for different retardation 
parameters e, analogous to the van der Waals liquid-gas isotherms. The extremities 
of the stress plateau, for different e, form the binodal curve. Above e = 1/8 (critical 
point) the constitutive curve is monotonic. The homogeneous branches and the 
plateau constitute the flow curve in planar shear, b) Limiting values of the stress 
allowing coexistence of bands, as functions of the retardation parameter e. a^^^ is 
the selected stress value at the interface. The values are consistent with the results 
in [21] obtained with a different method. 

The stress is expanded about the interface position, a = a* 

where a* = a{f*) is the value of the stress at the interface. By neglecting terms 

of order D^^"^ and higher we obtain, for all studied geometries: 

dlS,n = Sin - ^^^(1 - W,n) (3.3a) 

dlW.n = W.n - ~ ^^"^ (3.3b) 

Zin = (3.3c) 



For pipe flow we may also use X^n — Win, that together with Eq. (3.3c) implies 

We refer here only to single interface solutions. More complex solutions with 
an arbitrary number of interfaces can be treated in the same way. There are 
two types of single interface solutions, differing by the sequence of bands. Let 
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us refer to the solution with the high shear rate band at the left of the interface 
(towards smaller r) as (H — ), and to the solution with the high shear rate band 
at the right of the interface (towards bigger f) as ( — h). 

The inner and outer solutions should match at the interface (Prandtl matching 
principle, Ref. [34]), leading for the sequence (H — ) to: 

5+-(-oo) = sun W,t-{-oo) = (3.4a) 

^+-(00) = sun W,t-{oo) = (3.4b) 



or for the sequence ( — h) to: 



Sin 



Souti^ ] 



4-/(00) = sun 



W- 



-00 



w. 



out 



a 



wu{n = wu^* 



(3.4c) 
(3.4d) 



3.3 Equations of flow curves 



The flow curve expresses the relation between measurable quantities: velocity 
of the inner cylinder V and torque T in the Couette geometry, maximum 
velocity V and pressure gradient / in Poiseuille flow, gap velocity V and 
shear stress a in planar shear. 

In order to obtain the flow curves we use the shear rate conditions, Eqs. (2.8), 
(2.11), (2.13), and the relations (2.6b), (2.9), (2.12) between the total stress a 
and the position f. In the banded regime, the shape of the flow curve depends 
on the number of interfaces and on the sequence of bands. We shall show in 
the next section that in the presence of stress diffusion steady flow always 
corresponds or can be reduced by symmetry to a single interface. 

In Poiseuille flow, using Eqs. (2.6b), (2.8) or (2.12), (2.13), a single interface 
flow curve with the high shear rate band near the wall (right of the interface) 
is described by: 



y(/, n = 



/ 



7 {a)da + / 7 {a)da 

J a* 



(3.5) 



where a* — a{f*) is the value of the total stress at the position of the interface. 
The inverse sequence of bands is prohibited in Poiseuille flow, because the high 
shear rate band can not exist at f = 0, where a — 0. 

In Couette flow, using Eqs. (2.9), (2.11) a single interface flow curve with the 
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high shear rate band near the inner cyhnder (left of the interface) is given by: 



ir/(l+p)2 (7 Jct* (7 



/r/(i+p) 



(3.6) 



The flow curve for the inverse sequence of bands (high shear rate band near 
the outer cylinder) is given by: 



r/(i+p)2 



da 

7 (a)- 



^ ^ - , .da 
7 (^)^ 



(3.7) 



Planar flow represents a special case because the total stress is constant 
throughout the gap (Eq. (2.6a)) and thus: 



V{u, a) = wy (a) + (1 - i^)7 (o") 



(3.8) 



where is the proportion of high shear rate band in the flow. As we shall 
argue below, the flow curves above are well defined because the value it* of 
the shear stress at the interface is a geometry- independent constant. 



4 Uniqueness and stability of the steady flow 



4-1 An existence conjecture for the inner layer solution 



The inner layer solution controls the existence and stability of steady banded 
fiows. Because we are dealing with singular perturbations, this control can 
be performed even by a very thin interface (very small D). Before stating an 
important property of the inner layer solution, we note that steady banded 
flow is a particular case of a moving interface solution f*{t). We may look for 
matched asymptotic solutions in this case as well and one may show by using 
the change of variable 

f^t^ (4.1) 



that the moving interface inner layer solution should obey: 

d^.Sin + cdfSin = Sin - ^-^^(1 - W^n) (4.2a) 
d'.mn + cdrWin = - '^^"^'^^ ~ '^^"^ (4.2b) 
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where 

(4.3) 

is the rescaled velocity of the interface. 

Conjecture. There is a unique value a* = dgei G {o'tottomi^) , o'topi^)) of the 
total shear stress at the position of the interface, such that a steady inner layer 
solution for the banded flow of the JS-d model (Eqs. (3.3) exists obeying the 
Prandtl matching principle, Eqs. (3.4). Moreover, moving interface inner layer 
solutions obeying the Prandtl matching principle for Eqs. (4-2) exist for val- 
ues of a* in a neighbourhood of dgei, with velocities c{a*) that are well defined 
functions of a* . The stationary solution satisfies c{&sei) — and, further- 
more, > when the high shear rate band is at the left (towards smaller 
values off) of the interface, and -^la^^i < for the inverse sequence of bands. 



Remark 1. The conjecture considers the inner layer solution, which obeys 
the same equation for all flow geometries. This implies that the selected 
value of the stress is geometry-independent. After the rescalings, the dif- 
fusion coefficient does not enter Eqs. (4.2). Thus, the selected stress is in- 
dependent of V. As shown by Eq. (4.3) the speed of the interface df*/dt 
scales with i)^/^, thus the role of D is to speed up kinetics of non-steady flow. 

Remark 2. A stationary inner layer solution which obeys the Prandtl match- 
ing principle represents a heteroclinic orbit of the 4D dynamical system 
l^, c^fS", H^, connecting the hyperbolic flxed points ^S~^{asei), 0, 

W~^{asei), o| and ^S~{asei),0,W~{d-sei),Oy As discussed in [19,21], and as 
is valid rather generally in these cases, the value of the parameter a* — &sei 
allowing this solution is isolated (actually unique in the interval {(Tbottomi^), 
atopi^)))- The condition on the sign of the derivative means that an 
increase of the stress at the interface above a^ci produces a movement of 
the interface that decreases the size of the low shear rate band. Conversely, 
a decrease of the stress at the interface below agei produces a displacement 
of the interface that increases the size of the low shear rate band. The 
same property can be found for other reaction-diffusion systems {e.g. the 
FitzHugh-Nagumo model of nerve conduction in biophysics) and has been 
occasionally referred to as "dominance principle" [35,36]. 

Remark 3. We have numerically checked this conjecture for various values 
of the parameter e of the JS-d model (see next Section). All of these features 
are easy to prove for toy models that lead to integrable dynamical systems 
[22]. 
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4-2 Numerical test of the conjecture 



In order to test the conjecture for the JS-d model, and determine the relation 
between c and a* for different values of the unique parameter e of the model 
in reduced variables, we have numerically integrated the following system of 
partial differential equations: 

dtS = l3dlS + c^dfS -S+ ^^^^ ~ ^ {l-W) (4.4a) 
dtW = l3dlW + c^dfW -W+ '^^^^^^ ~ (4.4b) 

with 



SiL) = S-ia,oUoM), Si-L) = S+iatopie)) (4.4c) 

W{L) = W-{a,oUom{^)), W{-L) = W+{atop{e)), (4.4d) 

where 

r + L 

cr{r) = (7top(e) + ((J6ottom(e) - crtop{e))^j^ (4.5) 

is conveniently chosen to scan the interval [o'bottom{^),c!'top{^)]- Provided that 
the length L is taken much larger than the interface width, i.e. \^ <^ L, the 
stationary solution is an interface in a position corresponding to a value of the 
stress equal to a* = (3"*(c). j3 can be interpreted as the diffusion coefficient and 
cr(f) as the total stress inside the gap of width 2L although the applicability of 
the method does not rely on its connection to a concrete physical problem, but 
on its formal analogy to Eqs. (2.5) and (2.7), in the asymptotic {L/^/P — > oo) 
limit. 

In this way we determine aseii^) = cr*(c = 0; e) and -^{(Jseh e) = (^) (c ~ 
0; e). The first dependence will be used in calculating the fiow curves in Section 
5, while the sign of the latter constitutes a proof of the dominance principle. 

In order to determine precisely the position of the stationary interface and 
hence the selected stress, we have solved Eqs. (4.4) for several values of (3. 
Changing (3 modifies the interface width, analogous to a change of the diffusion 
coefficient, but does not alter the selected position of the interface for small j3. 
The selected position of the interface and the corresponding value of the stress 
are given by the cross-over of several profiles obtained for different values of 
(3 as in Figs. 2a,b. 

This algorithm is more convenient than the one used in [21] that considers a 
constant stress a throughout the gap (planar shear) and tunes "by hand" the 
value of a to obtain a stationary interface. Here the tuning of the stress is 
automatically obtained while integrating Eqs. (4.4). 
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Fig. 2. a) Strain rate 7 as a function of position for three different stationary so- 
lutions (c = 0) to Eq. (4.4) for /3 = 2, 4, 8. (b) Stress as a function of position. 
The cross-over of several stationary profiles gives the selected stress at the position 
noted by the vertical line, c) Derivative at = of c"' [a* ) as a function of the 
retardation parameter e, for an interface with the high shear rate band at its left. 
Here, cTsei = 0.4799, compared to cTseZ = 0.4827 found in [21] for the JS model in 
planar Couette flow by a different method. 

The selected value of the stress is compared with abottom and atop in Fig. lb 
and agrees with the result of the method in [21]. 

In order to test the dominance principle we have calculated Changing 
f ^ — r and c — >• — c leaves Eqs. (4.4) invariant. Thus if c^~{a*) and c~^{a*) 
are the velocities of the interface with the high shear rate at its left and right, 
respectively, then c~'^{a*) — —c^''{a*). Hence, our check that > 0, as 
shown in Fig. 2c, represents the sufficient test of the dominance principle. 



4-3 Number of bands and stability of banded flow 



The numerically tested conjecture implies that the interface is stationary in 
a position inside the gap where the value of the stress is asei- This position 
is unique (f^ei = (f /o-se/)^/^ for Couette flow, f^e; = ^sei/ f for pipe flow) 
because in steady flow the total stress is monotonically dependent on the 
radius. Therefore, in Couette and pipe geometries banded flow has only one 
interface separating two bands. For Poiseuille sht flow, the dependence of the 
total stress on position is symmetric with respect to the middle of the slit, so 
there are two interfaces at fgei — '^^sei/ f ■ 

For both Poiseuille slit and pipe flow the sequence of bands is unique, with 
a centre low shear rate band and outer high shear rate bands (or an annular 
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band in pipe flow) . In Couette flow, two sequences of bands are possible under 
certain conditions that we shall discuss in the next section. There are no 
restrictions on the sequence and order of the bands in planar flow because the 
total stress is constant throughout the gap; interfaces can exist anywhere, and 
any number and sequence of bands is possible, provided that the bands are 
much thicker than the interface width. 

The stability of banded steady flow in presence of diffusion terms was tested 
in the Couette geometry by numerically evolving the dynamical equations [19] 
for both normal (high shear rate at the inner cylinder) and inverted sequences 
of bands. We have also proved (Appendix A) the stability of different types of 
banded solutions in Poiseuille and Couette ffows with respect to perturbations 
of the position of the interface. A linear stability analysis for ID perturbations 
around a banded proflle, as well as the stability with respect to surface waves in 
2D, involves lengthy functional analysis and will not be discussed here. Linear 
stability of banded flow for ID perturbations in slit Poiseuille geometry was 
shown in Ref. [24] for a simplified model with only one order parameter. For 
the JS model (without diffusion) linear stability of the banded flow with any 
flnite number of bands was proven in [15] for ID perturbations in Poiseuille 
geometry, and possible instabilities induced by surface waves were discussed 
in [11]. 

Although linearly stable, the inverted sequence of bands is metastable because 
the interfaces separating a nucleus of low shear rate band within the high shear 
rate band will flnd themselves in a region of stress lower than &sei and according 
to the dominance principle will have non-zero velocities oriented such that the 
nucleus will grow [22] . The same is true for a nucleus of high shear rate band 
in the region where the total stress is higher than dsei- 



5 Applications of the results 

5. 1 Flow curves and eoctremities of branches 

Flow curves can be obtained from Eqs. (3.5-3.8). These curves have several 
branches corresponding to homogeneous and banded flow. The extremities of 
different branches are flxed by the following conditions: 

(1) Homogeneous low shear rate flow is possible if 

0<a{f) <atop{e),W (5.1) 
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(2) Homogeneous high shear rate flow is possible if 

cr{r) > hottom{e),'ir 

(3) Banded flow is possible if 

for some position r*^i inside the flow. 

The sequence (H — ) is allowed if, in addition: 

< a{f) < atopie),\/f, r > f^^i 
cr{r) > hottom{e)yf, f < f*^i 

while the sequence ( — h) is allowed if, in addition: 

< a{f) < (7top(e),Vf, f < rl^i 
o-(r) > <3-6ottom(e),Vf, f > f*^i 



(5.2) 



(5.3) 



(5.4) 
(5.5) 



(5.6) 
(5.7) 




Fig. 3. Shear rate profiles (thick line) for banded flow (to each value of the total 
stress a corresponds a position r inside the gap) in a) Poiscuillc geometry, b) Couette 
geometry p = 0.05, e = 1/20 (region I of Fig. 5), with two possible sequences of 
bands ((H — ) solid, high shear rate band at the inner cylinder), and ( — h) dashed, 
high shear rate band at the outer cylinder), c) Couette geometry p = 0.1,e = l/20 
(region IV of Fig. 5 inverted sequence is not possible). The interface was considered 
of negligible width, hence it is represented by the segment a = (Tgei- 

In Poiseuillc and Couette geometries, any banded steady flow profile can be 
conveniently represented as a segment on the local constitutive curve 7(0"), a 
solution of Eqs. (3.1) for variable stress a (Fig. 3). Using the relation between 
f and a (Eqs. 2. 6b, 2. 9, 2.12) one can map a steady state 7(f) onto a segment 
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of the local constitutive curve 7(0") for stresses < a < f (Poiseuille), or 
r/(l + < a < r (Couette). The outer layer solutions lies on the local 
constitutive curve, while the inner layer solution is a horizontal segment at 
a — a sell corresponding to a small width interface. 

5.1.1 Planar shear 

The equation of the low shear rate branch is y = 7 (a), with < a < (5"top(e), 

and the high shear branch can be described as V = 7^(0"), with (Jbottom,{^) < 
a. The banded flow corresponds to the plateau a — Ogei- The homogeneous 
branches of the flow curve are on the local constitutive curve represented in 
Fig. 1. 

5.1.2 Poiseuille flow 

The low shear rate branch can be obtained from Eq. (3.5) by choosing a* = f 
(interface near the wall at f* = 1), and according to Eq. (5.1) with the restric- 
tion < / < atop{e)- The unique banded profile ( — h) (low shear rate at f = 0) 
is described by Eq. (3.5), choosing a* — dgeu according to Eq. (5.3) with the 
restriction / = asei{^) /f*^i > o"se/(e)- Then Eq. (5.7) is automatically fulfilled 
because it becomes equivalent to (Xhottomi^) < (^seii^) < o'topi^), that is true for 
e < 1/8 (actually for all values of e for which the local constitutive curve is 
nomonotonic. Fig. 1). Eq. (5.5) is not fulfilled (thus the inverted sequence of 
bands is forbidden) because it becomes equivalent to fr > Cbouorm^'f^ < "^tei^ 
obviously not true for f = 0. The banded branch continues until / = 00 and 
f* = 0. The low shear rate band is continuously squeezed until it disappears 
in the limit / = 00. Thus, the flow curve has no high shear rate branch, be- 
cause the high shear rate band occupies the entire volume only at / = 00. The 
flow curve in this geometry is represented in Fig. 4a for different retardation 
parameters e. 

5.1.3 Cylindrical Couette flow 

In Couette flow Eq. (5.5) is always fulfilled and the normal sequence ( — h) 
of bands is allowed for all shear rates along the plateau (see Fig. 3b-d). Nev- 
ertheless, the inverted sequence ( — h) is not allowed for the values of T,p,e 
when the total stress at the inner cylinder is higher than atop, or when the 
total stress at the outer cylinder is lower than abouom, see Fig. 3c. A detailed 
analysis of this condition given in Appendix B shows the existence of five 
possibilities and regions in the plane of parameters (p, e), explained in Fig. 5. 

Some examples of flow curves in the Couette geometry are presented in Fig. 4. 
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0.00 



0.55 



0.45 



0.35 




Fig. 4. Flow curves in reduced variables for a) Poiseuille flow, b-d) Couette geometry 
with p = 0.01 (b), p = 0.06 (c), and p = 0.11 (d). For Couette flow, 1/+" and 
V~'^ correspond to flow curves with the high strain rate band near the inner and 
outer cylinders, respectively, (b) and (c) correspond to Regions I and II, while (d) 
corresponds to Region IV in Fig. 5. 



Eq. (A. 8a) implies that the branch (h — ) has a positive slope |£r > 0, while 

Eq. (A. 8b) implies that the branch ( — h) has a negative slope |^ < 0. The 
latter does not imply a mechanical instability, because any mechanical drift 
is constrained by the imposed mean shear rate. Furthermore, the stability of 
this sequence of bands is shown by our numerical simulations [19] and by the 
stability analysis with respect to the position of the interface (Appendix A). 
As suggested in [16,17], stability of banded flow is not directly connected to 
the slope of the flow curve, but rather to the slope of the local constitutive 
curve of the homogeneous bands; if the latter is negative for at least one band 
(intermediate branch in Fig. la), then the flow is unstable. 

Fig. 4b corresponds to values of the parameters (p, e) chosen in Region I of 
Fig. 5, with both sequences of bands allowed. Fig. 4c corresponds to param- 
eters chosen in Region II (the inverted band sequence is forbidden at higher 
strain rates), while Fig. 4d corresponds to parameters chosen in Region IV 
(shear rates allowing the existence of the inverted sequence are limited both 
at upper and lower values). Interestingly, the branch corresponding to the in- 
verted band sequence no longer crosses the normal band sequence branch for 
p = 0.11, e = 0.05 in Fig. 4d. Note that the "plateau" in the banded region is 
steeper (larger dV /dV) for a more highly curved geometry (larger p). 
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Fig. 5. Existence conditions for the inverted sequence of bands in Couette flow, 
separating 5 regions I-V in the plane {e,p). The dots • correspond to the flow 
curves in Fig. 4. Inverted bands arc only possible in Regions I (for all shear rates 
along the plateau), II (for a domain of shear rates truncated at the right end of the 
plateau) , III (for a domain of shear rates truncated at the left end of the plateau) , 
and IV (for a domain of shear rates truncated at both ends of the plateau). 

The various dangling segments corresponding to different branches of the flow 
curves may be the origin of hysteresis phenomena. Part of these phenomena 
were described in a previous paper [19] using numerical integration of the dy- 
namical equations. We have shown there that some segments of the flow curve 
can be reached only by special experimental scenarios. This is in agreement 
with the results of this work. For instance, starting from rest in Fig. 4b, when 
e = 1/25, the system can flrst reach the branch (H — ) within the interval of gap 
velocities Vi < V < V2 and the branch ( — h) within the interval V2 <V <V3 
(Fig. 4b). In order to reach other parts of the banded branches one must, for 
example, prepare the system with a given band sequence and then adiabati- 
cally change the value of the gap velocity in order to scan the entire length of 
the branch. This is possible for p < 0.1 for both bands (H — ), but for p — 0.11 
in Fig. 4d all start-up preparations of the flow end with the sequence (-1 — ) and 
it is impossible to reach the isolated branch ( — h). Numerically one could start 
from rest with a smaller value of p, reaching thus the ( — h) branch, and then 
change p adiabatically, but this is not a conceivable experimental procedure. 
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5.2 Interface Width 



Linearising Eqs. (3.3) at r — ±00 we obtain 

di\^^r]=M^i'-r I (5.8) 



1 SSin 1 














[5WJ 



where 5 Sin, SWin are small deviations of stresses with respect to asymptotic 
steady values and 

, ,l-[l-W^{^sei)]/e[a,,i-S^{asei)]/e\ 
M {asei) = I . I (5-9) 



Asymptotically, the interface profile can be approximated by a combination 
of exponentials, and the widths of the interface towards the high shear rate 
band and towards the low shear rate band are the characteristic decay 
lengths of the slowest of the exponentials on the two sides of the interface, 
related to the eigenvalues of M='=: 

\/Vt 

= r /r— (5-10) 

min|Re \^^Jxi (<^*,e)J| 

where i = 1,2 correspond to the two eigenvalues Xi,2 of M~ and to the two 
eigenvalues Xi,2 of (see Fig. 6a). The dependence of the interface width on 
the retardation parameter e is shown in Fig. 6b. The interface is asymmetric 
(thicker in the low shear rate band). 



The non- vanishing imaginary parts of \jxi,2 imply damped oscillations of the 
interface profile towards the high shear rate band. The wavelength of these 
oscillations is the inverse of llm^x^l and because |lmy^x^| < Min(Rey^x^) 
(Fig. 6a) the width of the interface is smaller compared to this wavelength the 
overdamped oscillations are hardly noticeable in Fig. 2a showing the interface 
profile. 

As seen in Fig. 6a, a bifurcation occurs for low polymer viscosity. For e > 0.1, 
the eigenvalues Xi,2 no longer complex conjugate (Rex| ^ Kext)- This 
bifurcation produces a discontinuity of the first derivative of w^{e) at e = 0.1. 
At the same time the imaginary parts of X12 vanish. This particularity of the 
JS model was also noticed in [14]. At e = 1/8, which is the limit of existence 
of banded solutions (above e = 1/8 the local constitutive curve is monotonic) 
Xi 1 Xi vanish, therefore both widths diverge. Fig. 6b. 
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Fig. 6. (a) Eigenvalues of M^((Tsei) the matrix of the hnearised steady flow equa- 
tions, for the high (+) and low (-) shear rate bands and (b) interface width on 
the high or low shear rate side of the interface, as a function of the retardation 
parameter e. 

The interface width scales like \/Vt, which is typically of the same order as 
the polymer chain size [22] . The prefactor can be rather big for small polymer 
viscosity (it diverges at the critical point e = 1/8 like (1/8 — e)~^/^) and it is 
geometry independent. Thus, the width of a thin, steady interface should be 
the same in planar shear, as in Poiseuille and Couette flows for any curvature. 
Ref. [30] determined the width of the interface for micelles experiments in 
cone and plate geometry, and reported an increase in width of a factor 3 when 
the cone-plate angle changes from 4" to 7°. Although we did not investigate 
this particular geometry here, we expect that rather generally the width of 
a thin, steady interface depends only on the constitutive model and does 
not depend on the geometry. The widths reported in [30] are much larger 
than the polymer chain size, therefore they are most probably produced by a 
different mechanism. A possible explanation is the dynamical broadening due 
to travelling surface waves in two-layer flows [37,38] . 



6 Conclusions 



We have discussed the steady banded flow of the Johnson-Segalman model, 
considering ID shear rate profiles in Poiseuille, cylindrical Couette and planar 
shear geometries, under controlled shear rate conditions. The introduction of 
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a non-local diffusion term in the JS local constitutive relation lifts the contin- 
uous degeneracy and the history dependence of the steady flow solution for 
Poiseuille and cylindrical Couette geometries, and generally for any geome- 
try that imposes a non-uniform total shear stress. In order to analyse these 
phenomena we used matched asymptotic techniques which are common in the 
field of reaction-diffusion systems describing processes like combustion or nerve 
conduction, but which, to our knowledge, are new in the field of creeping fiow 
of complex fiuids. The use of the JS model in this paper should be considered 
only as an example of application of these methods in order to investigate the 
consequences of diffusion terms on banded flow. 

We showed that in Poiseuille and cylindrical Couette geometries steady banded 
flow has two bands, which seems to be the case in birefringence experiments 
[31,32]. NMR visuahsation of more than two bands as reported in Ref. [29] in 
the Couette geometry could be due to slowly decaying (or pinned) transients. 
Note that the speed of a moving interface scales as D^^'^. Another possibil- 
ity, not taken into account here but that needs further investigation, is that 
concentration effects may stabilise banded flows with more than two bands. 
Thus, a depletion layer at the inner cylinder increases locally the selected 
stress and may equilibrate a second interface closer to the inner cylinder than 
the first equilibrium position which corresponds to higher concentration and 
lower selected stress. 

In the absence of diffusion steady banded flow can form any number of bands. 
The same arbitrariness of the number of bands was shown in steady Poiseuille 
flow [16,17]. In numerical simulations of Ref. [25] no selection was noticed, 
except for the one imposed by the mesh size. Recently, we showed numerically 
[19] that even in Couette flow, solving the flow in the absence of diffusion is 
an ill-posed problem, unstable to noise in the initial conditions. The lack of 
robustness of the flow increases when the curvature decreases, in the sense that 
larger and larger regions inside the gap may split into multiple bands under 
perturbations of the initial conditions. Refs. [12,39] found an unique, single 
interface solution even in the planar shear of the JS model without a diffusion 
term, although this interface has a flnite width that suggests diffusion terms 
intrinsic to the finite element scheme. With diffusion the system of dynamical 
equations becomes parabolic leading to compactness of the global attractors 
and to stability with respect to initial conditions [24]. 

In certain conditions in cylindrical Couette ffow (low curvature and high poly- 
mer viscosity), both normal and inverted band sequences are allowed for the 
same shear rate. Nevertheless, the inverted band is metastable and can be 
reached from rest only on the second half of the plateau [19]. Although this 
branch has not been observed experimentally, its existence and linear stabil- 
ity is proven rigorously and may have consequences on the kinetics of shear- 
banding within a restricted domain of control parameters. 
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An interesting feature shown by experimental flow curves and reproduced by 
our theoretical model is the presence of a selected stress plateau, reminiscent 
of isotherms at liquid-gas or decompositions of binary systems first order phase 
transitions. In thermodynamic transitions the plateau is fixed by a Maxwell 
type construction, justified by the equality of the chemical potentials of the 
coexisting phases. In the model we present here the equilibrium of "phases" 
is a condition of existence of a stationary interface, and relies on the balance 
between stress diffusion across the interface and non-linear stress relaxation. 
The retardation parameter plays the role of temperature, and its critical value 
e = 1/8 is analogous to an equilibrium critical point. The extremities of the 
plateau for varying e describe a curve that is analogous to a binodal for phase 
separation, the critical point being the top of the binodal (Fig. la). As in first 
order transitions the width of the interface separating the two bands diverges 
at the critical point. Whether the studied phenomena have more than a formal 
connection to thermodynamics is still an open question. 

Apart from settling some basic questions concerning how the presence of dif- 
fusion terms affects measurable flow curves, we hope that the results pre- 
sented here will stimulate more accurate and diversified experiments, eventu- 
ally searching for various unusual features such as inverted band sequences 
and the special kinetic pathways involved. This could be a test of the theo- 
retical assumptions but could also increase the operational range in view of 
possible applications. 
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A Stability of the flow with respect to changes in the position of 
the interface 

Stability with respect to the position of the interface means that any displace- 
ment 5r* of the interface leads to a moving interface solution with velocity 
oriented towards the equilibrium position fseu i-^-'- 

{dr*f < (A.l) 

V 



dc 

c [a {fsei + Sr*)] 5f* = — 



<T„„I 



dd_ 
dr* 
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The partial derivative occurring in Eq. (A.l) can be obtained as 



da* 



df* 



V 



dV 




df* 


a* 


av 




aa* 


f* 



Using Eqs. (A. 2, 3.5, 2.6b) for Poiseuille flow we obtain: 



dv_ 

df 



f 



da* 


av 

_ df 


5-* af* 


a* _ f 


df* 


av 

V Q^* 


I av 

r 9f 


a* 9a* 


1-P{f,asei) 

f* 



(A.2) 



(A.3) 
(A.4) 



where 



sell 



f[l^{f)-n^sel,f)] 



(A.5) 



According to the tested conjecture, an increase of the total stress will displace 
the interface towards the middle of the slit or pipe < 0), so using Eq. (A.4) 
the stability condition (A.l) for plug flow in Poiseuille geometry reads: 



P{f,asei)<l. 



(A.6) 



Using the fact that 7 (a) is an increasing function of a it follows from Eq. (3.5) 

that V{asei, ^sei) < 7 {c^sei) < 7 {c^sei) and therefore P{aseu c^sei) < 1- P{f, ^sei) 
is a monotonically decreasing function of /, provided that 7 {a) is increasing 
with a, as shown by 



dP 
df 



asel 



d^\f) 

df 



1 {f)-V{aseuf) 



< 0. 



(A.7) 



Because the lowest allowed value of / for banded flow is f — asei, the stability 
criterion is always satisfled. 



A similar stability analysis can be performed in the case of the cylindrical 
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Couette flow. Using Eqs. (A.2), (3.6), (3.7), (2.9) for Couette flow we obtain: 



dV+- 



dr 

dV-+ 



or 



da* 



d-r" 
da* 



df* 



v+- 



V- 



7 (r)-r(r/(i+p)^) 
2r 

-^^(r/(i+p)^)+r(r) 

2T 

l-C+-{t,asei,p) 
l-C-+(f,a,e/,p)' 



(A.8a) 
(A.8b) 
(A.8c) 
(A.Sd) 



with 



C+ {t,asei,p) 

C~+{t,asehP) = 



7 i^sel) - 7 (g"ge;) 

"^+(f)-r(f/(i+p)2) 

7 (g"^e;) - 7 i^sel) 
^+(f/(l+p)2)-f (f)' 



(A.8e) 
(A.Sf) 



where the sequence H — refers to the high shear rate band near the inner 
cylinder and the sequence — h refers to the high shear rate band near the 
outer cylinder. 



The stability criteria for the two sequences are respectively: 

C+-(f,a,e/,p) < 1 
C-+(f,a,e/,p) > 1 



(A.9a) 
(A.9b) 



Because f/(l +p)^ < asei < T, and provided that both 7 (a), and 7 (a) 
are increasing functions of the total stress a, the stability criteria are always 
fulfilled for both sequences of bands. 

To conclude, banded flow is stable with respect to displacements of the in- 
terface in Poiseuille sht and pipe geometry. In cylindrical Couette flow, both 
sequences of bands are stable with respect to displacements of the interface, 
provided that their existence is allowed by further criteria discussed in Ap- 
pendix B. 
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B Existence conditions for band sequences in cylindrical Couette 
geometry 

The homogeneous low shear rate branch is obtained from Eq. (3.6) with a* = F 
(interface at the inner wall f* = 1 ), and according to Eq. (5.1) with the 
restriction < T/P < atop{^),^r G (1, 1 +p), which is equivalent to < F < 
(3"top(e). The high shear rate branch can be obtained from Eq. (3.6) with a* = 
r/ (1 +p)^ (interface at the outer wall r* — 1+p), and according to Eq. (5.2) 
with the restriction f/f^ > atop{e),yf e (1, 1+p), i.e f > (1 + py&bottomie)- 

Banded flow with the two possible band sequences is described by Eqs. (3.6, 
3.7) with a* = dgei- Eq- (5.3) imposes the limiting values of the torque: 

Usel < f = &sel{r*sei? < ^selO-+pf (B.l) 

The two limiting values correspond to the steady interface touching the inner 
(f*g; = 1) and the outer (f*g; = 1 + p) cylinder, respectively. 

If Eq. (5.3) is fulfilled, Eq. (5.5) is automatically fulfiled because it becomes 
equivalent to crbottomi^) < (^sei{^) < o"top(e)- Thus the sequence of bands (H — ) 
(high shear rate band at the inner cylinder) exists for all torque values such 
that: 

tel^:^{asei:asei{l+pf). (B.2) 
If Eq. (5.3) is fulfilled, the condition 5.7 becomes equivalent to: 

f e /r := ((1 + pfabottom: ^top): (B.3) 

hence the ( — h) sequence of bands (high shear rate band at the outer cylinder) 
exists for torques F e = 7p fl /p- The interval /_+ can be either empty, 
or strictly included in Jp, or coincide with /p. There are several regions in the 
plane of parameters (e,p), corresponding to these situations. 

(I) 0-6ottom(e)(l+p)^ < ^seh ^top{^) > ^seli'^+pf: ^bottom{^){'i- + pf < ^top{^)- 

In this region both band sequences are allowed for torques in the full 
interval Jp because Jp C It and therefore / ^ = /p. 

(II) (3-6ottom(e)(l + pf > CTsel, Crtop{e) > (Jselil + ^»)^ 0'feottom(e)(l +p)^ < 

(Jtop(e). In this region the sequence (H — ) is allowed for torques within Jp, 
while the sequence ( — h) is allowed in the interval /_+ — {d'bottom{^){^ + 
p)^, <5"sei(^)(l obtained by truncating /p at its low stress limit. 

(Ill) (T6ottom(e)(l+p)^ < ^sel, ^top{^) < CTselil+pY, (^tottom (e) (1+ P)^ < ^top{e) . 

In this region the sequence (H — ) is allowed for torques within Jp, while the 
sequence ( — h) is allowed in the interval /_+ = {d'sei{€-),d'top{€)) obtained 
by truncating /p at its high stress limit. 
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(IV) abottom{e)(l+py > asei, ^top{e) < asei{l+pf, <3-6ottom(e)(l+p)^ < o-top(e). 
In this region the sequence (H — ) is allowed for torques within Jp, while the 
sequence ( — h) is allowed in the interval /_+ = (o"ftottom(e)(l ^top{^)) 
is obtained by truncating /p at both limits. 

(V) &bottom{e){l+pf > ^seh ^top(e) < Crsel{l+pf-, Crbottom{e){l+pf > ^top{^)- 

In this region the sequence (H — ) is allowed for torques within /p, while the 
sequence ( — h) is forbidden for all values of the torque, because /_+ = 0. 

These five regions are shown in Fig. 5. Region III is very narrow, and prac- 
tically represents the frontier between the Regions I and IV. In general, for 
given r and e the inverted band sequence is forbidden for sufficiently highly 
curved geometries (large enough p). A less formal approach to obtain the same 
results is shown in Fig. 3. Fig. 3c corresponds to the case of large curvature 
(region IV), in which the value of the stress does not allow for the existence 
of the low shear rate band at the outer cylinder forbidding thus the sequence 
(-+), because r/(l + p)^ < (Tbottom- For small curvature (Fig. 3b) the values of 
the stress at the inner and outer cylinders allow the existence of both high and 
low shear rate bands, thus both sequences of bands can exist. The possibility 
of an inverted band sequence was also found in Ref. [40] for a two fluid model 
of shear banding in polymer blends. 
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